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The condensation energy of the homogeneous electron gas is calculated within the density func- 
tional theory for superconductors. Purely electronic considerations include the exchange energy 
exactly and the correlation energy on a level of the random phase approximation. The singlet su- 
perconductivity is assumed, and the Coulomb interaction is studied with a model pairing potential 
at the angular momentum up to 1=9 and at densities l<rs<10. The homogeneous gas remains 
non-superconducting up to rs~9. Very weak negative value of the condensation energy has been 
found for f-waves and higher-^ pairing at rs = 10. 
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I. INTRODUCTION 

For superconductivity both the electronic and the 
electron-phonon interactions are important. While for 
the BCS 'li| superconductors the electron-phonon con- 
tributions dominate, for the novel high-Tc materials the 
Coulomb interactions play a key role. The Eliashberg 
theory 01 , being a generalization of Migdal's theorem Q , 
is able to treat the electron-phonon interaction in both 
the weak and the strong coupling regime, but the elec- 
tronic interactions are averaged to a single parameter, /u* . 
This might be insufficient even for the strongly correlated 
systems 4]. In this work, we focus on purely electronic 
interactions in the superconducting homogeneous gas for 
which the results obtained in last forty years still remain 
controversial. 

The existence of superconductivity at higher angular 
momentum pairing without phononic contributions was 
suggested by Kohn and Luttinger jSJ], in 1965. The mech- 
anism proposed there was based on the presence of the 
long-range oscillatory potential in ordinary space due to 
the sharpness of the Fermi surface, and the fact that 
Cooper pairs Q could form taking the advantage of the 
attractive regions. Interestingly, some features of the 
phonon spectra have been explained due to Friedel os- 
cillations 0]- 

More than twenty years ago, Takada p solved the 
Eliashberg equations and estimated the transition tem- 
perature, Tc, due to the plasmon exchange. His solu- 
tion assumed weak electron-phonon coupling for which 
the Kirzhnits, Maksimov, Khomskii approximation [9j 
can be justified. Other authors, Rietschel and Sham 
|lO] and Shuh and Sham [ij , solved the strong coupling 
Eliashberg equations linearized in the gap function. For 
the Coulomb interactions, they also assumed the ran- 
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dom phase approximation (RPA) [l3| and found unre- 
alistically high critical temperatures. Later papers in- 
cluded vertex corrections yielding a counterpart for the 
plasmon exchange which overestimate the Coulomb at- 
traction. Obtained this way, critical temperatures were 
much lower [H 111 El El . 

Nowadays, we can use a standard solid state method, 
the density functional theory [13 (DFT), generahzed to 
the superconducting state (SCDFT) by Oliveira, Gross 
and Kohn [l3| in 1988. The formal framework, which 
we briefly describe in the next section, has been devel- 
oped over more than ten vears Il9t and extended for 
the relativistic superconductors ^(j\. The first solutions 
of the SCDFT scheme for simple metals predicted the 
critical temperatures |23 quite well, although the elec- 
tronic correlations have been taken into account only by 
the Thomas-Fermi screening of the Coulomb interaction. 
Then, the semi-phenomenological correlation functional 
was constructed and tested, firstly for Nb l22ll. and later 
for YBCO 23] using the eight-band model Hj. The for- 
mulation of the local density approximation (LDA) for 
superconductors was given by Kurth et al. |25j . in 1999. 
The attempt to construct the exchange-correlation func- 
tional from first principles was made |25j using as a start- 
ing point the RPA 12J generalized to the superconduct- 
ing state |2y| . The condensation energy of the homoge- 
neous gas has been calculated within that scheme for the 
model pairing potential of s-type |23| and no supercon- 
ductivity has been found up to rs=5. 

After the discovery of the anisotropic gaps in B2212 
[271 1 and YBCO |23| by the angular-resolved photo- 
emission experiments, it is interesting to look closer 
at the higher angular momentum channels within the 
SCDFT method. This is a step toward DFT calcu- 
lations for the superconducting state of high-Tc com- 
pounds. Recently, the condensation energy for systems 
with the anisotropic gaps has been studied by Haslinger 
and Chubukov within the Eliashberg theory adapted for 
the spin-fluctuation mediated pairing 29, 30J (see also 
references therein). 



The fact that, the model calculations based on re- 
sults obtained for a homogeneous gas have been per- 
formed for strongly-correlated layered superconductors 
by Kresin with co-workers [31| and Seibold [S^, shows 
that our studies can contribute to understanding physics 
of complicated systems. The data for a homogeneous 
gas obtained with our code [SJl can also give basis for 
the parametrization of an LDA functional for inhomo- 
geneous superconductors in a simil ar p hilosophy as an 
ordinary LDA is a parametrization |34| of the quantum 
Monte Carlo data J33 calculated for a representative set 
of Ts (Wigner radius) values. 

This paper is organized as follows: after a brief descrip- 
tion of the theoretical background in the next section, we 
show the results of the condensation energy calculations 
in section III, and make a comparison with findings of 
other authors in section IV. 



II. THEORETICAL BACKGROUND 

The framework of a DFT for superconductors was for- 
mulated by Oliveira, Gross and Kohn ilSj. It rests on 
1-1 mapping between the density 



and the superconducting order parameter 



(1) 



(2) 



on one hand and the electrostatic and pairing potentials, 
Ws(r) and As(r, r'), on the other hand. Here, we assumed 
the singlet pairing, however the triplet pairing can be 
treated analogously 36j. The noninteracting Kohn-Sham 
potentials are functionals of both the normal density and 
the order parameter: 



Vs[n,x]{r) = wo(r) + 



^3^/ '^(r') 



A,Kx](r,r') = Aoir,r') + ld'i 

A^c[n,x]i'r,r'), 



|r-r'| 



(3) 



(4) 



where wo(r) is the lattice potential, and Ao(r, r') is an ex- 
ternal pairing potential of an adjacent superconductor. 
The second term in As(r, r') is the anomalous Hartree 
potential. The exchange-correlation potentials are de- 
fined as functional derivatives of the xc free energy func- 
tional -F'xc["', x] over the normal and the anomalous den- 
sity: 



Vxc[n,x\iT^) 



A,eKx](r,r') = - 



SF^c[n,x\ 
Sn{r) 

SF^c[n,x] 



Sx*ir,r') 



(5) 
(6) 



The corresponding Kohn-Sham equations have the form 
of the Bogoliubov-de Gennes equations |18l l37| (/i is a 
chemical potential of the superconductor): 



Mfc(r) 



Wfc(r) 



+ ^;s(r) -/i 



2 

dV'A,(r,r')«fc(r'), 

—-+Vs{r)-n 



Ufc(r) 



Vk{r) 



d^r'A,(r,r')wfc(r'), 



(7) 



(8) 



and result from the diagonalization of the noninteracting 
Hamiltonian 



— + Vs{r)- fi 



(fr /dVA*(r,r')*T(i")*i(r') + -H^-c 



^a(r) 



(9) 



The exchange-correlation functional includes in gen- 
eral the electronic and phononic contributions 19]. Here 
however, we are interested in the electronic part only, 
treated within RPA [Tjl for the superconducting state 
|2g . The LDA scheme for superconductors J23 has been 
constructed analogously to the local spin density approxi- 
mation (LSDA). In superconductors, the order parameter 
plays a similar role to that of the spin-magnetization in 
LSDA. The electron gas is exposed to the external pair- 
ing potential of the superconductor, just as the LSDA 
gas is under the influence of a magnetic field. 

The exchange energy of the superconducting gas is 
given by the expression: 



f4^J.s,As] 
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1 tanh —Eu 
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(10) 



where E-^ ~ y/^t + |As(k)P is the quasi-particle spec- 
trum and ^k = 2"^ — ji. The anomalous Hartree 
energy, fAH , is a functional of the pairing potential only 



/Aff[Ais,As 



(fik (fik' 



in 



47 (27r)3(27r)3|k-k'|2 
As(k)A:(k') ^^ 



Ei^Ek, 
tanhr^EAtanh(^E)^,\. (11) 

The RPA energy results from the summation of bub- 
ble diagrams with normal and anomalous Green's 



functions |25L l2q| ■ G and F (see Fig.P), and can be written 
as 
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where the Matsubara frequencies, even i/„ — ^^ and 

odd ojn — a , enter the Fourier transform of the 
irreducible polarization propagator ns(q, ;/„), with q be- 
ing the momentum exchange of the interacting electrons, 
as follows: 



2 f d?k 



^[G(k,c^„)G(k + q,c^„ 



F(k,w„)Ft(k + q,w„ 



+ yn)]- (13) 



Evaluation of Green's functions in the polarization prop- 
agator leads to the explicit functional of the Coulomb 
and pairing potentials 
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We assumed above, that the Kohn-Sham orbitals are 
plane waves, because we are interested in the conden- 
sation energy of the homogeneous electron gas. The way 
to obtain the LDA functional for inhomogeneous super- 
conducting system is given in the Ref. |25|. Feynman 
diagrams for the contributions to the total energy con- 
sidered in this work are displayed in Fig.^and compared 
to those diagrams from earlier papers. 

The condensation energy can be obtained from 



Es-N = Tfei„[/Xs, As] — Tfci„[^s, 0] + fah[fJ-s, As 

/a;[^S,As] -/a;[/is,0] + 



(14) 



where the kinetic energy difference of the superconduct- 
ing and the normal state in a homogeneous gas at zero 



temperature is 



-^ kin [Ms : ^s 
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III. 



THE CONDENSATION ENERGY 
CALCULATIONS 



The first calculations of the condensation energy 
within the LDA for superconductors with the RPA func- 
tional were performed for the s-wave pairing only and 
no superconductivity was found 25]. We calculate the 
condensation energy of the homogeneous electron gas at 
zero temperature, assuming the model non-spherical gap 
function of the form 



Ar(k) 



5exp(^^-^)P™(k) 



(15) 



where 6 and a are parameters in units of /z and kp respec- 
tively, and P™(k) are associated Legendre polynomials. 
Above parametrization of the gap makes it possible to 
control the strength, the range and the angular shape of 
the pairing. In this work, we are mainly interested in the 
angular part. It will be clear from the further discussion 
that, the variational determination of S and a in such a 
way that the condensation energy is maximally negative 
would lead to either zero values for these parameters if 
Es^N is positive, or to infinite values of these parameters 
if Es-N is negative. As we will see later, the condensa- 
tion energy is also monotonic with I, but we were not 
able to predict this result from the analytical form of the 
expression (|14|l . 

The s-wave calculations appeared smooth in 6 
and /i parameters |25J in the range 0.01<cr<l and 
0.01<(5-100<1 at rs=0.1 and l<rs<5. Therefore in this 
work, we fixed the strength of the potential at (5=0. 01-/^ 
and the range of the pairing interaction at a—O.l-kp, 
and we present results for this choice of the parameters. 
Later, we will discuss changes in the condensation energy 
when it is calculated with two other sets of parameters 5 
and (7, namely with S—O.Ol-fi and a—O-Ob-kp, and with 
S^O.OOl-fi and a^O-l-kp- 

Turning to the details of the implementation, the eight- 
dimensional integrals of the energy functionals have been 
reduced by one dimension in electronic Matsubara fre- 
quency which can be evaluated analytically. Several 
singularities present in the formulas need special grids. 
For the calculation of /ah and f^ and the radial part 
of Jrpa, we used a modified Gauss-Legendre quadra- 
ture. For the angular part of Jrpa, we combined the 
Lobatto grid ^38] for the q-momentum integration, and 
the Sobolev's quasi random method [33 to generate the 
mesh used by the Monte-Carlo quadrature over the an- 
gular part of k-momentum. Details of the singularities 
and the parallel code are given in Ref. J33 . 



Kohn and Luttineer in 1965 
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Legend 




= bare Coulomb interaction 

— or » » = normal Green s function 

— or ^ » = anomalous Green s function 

— or < > or ^^v\^\* = summation 

= = self — consistent normal Green s function 



FIG. 1: Feynman diagrams for the total energy contributions considered by other authors IrLIhIkII and in this work. 



We will focus now on the results obtained with 
(5=0.01-/i and a—O.l-kp. We see in Fig. 2 the angular 
momentum dependence of the condensation energy and 
its components: the anomalous Hartree energy {/ah), 
and the diflFerence (S-N) of the exchange energy f^; be- 
tween the superconducting and the normal state and that 
difference (S-N) of the RPA correlation energy /rpa, cal- 
culated at rs = l and at rs=10. The condensation energy 
and all its components decrease monotonically with the 
angular momentum. The anomalous Hartree energy is 
the biggest positive component, almost completely bal- 
anced by the negative RPA energy which acts in favor of 
superconductivity. The exchange energy difference be- 
tween the superconducting and the normal state is pos- 
itive, and tends to destroy pairing. We show results for 
all momentum numbers I from to 9. However for the 
antisymmetric fermionic function only the even numbers 
make sense, because we assumed the singlet spin pairing 
for the order parameter. The condensation energies at 



rs = l are positive for 0<l<9. At the density rs = 10, the 
s-wave pairing also do not allow superconductivity. 

In order to show the energetics as a function of the 
density of a homogeneous electron gas, we present, in 
Table 1, the condensation energy and its components for 
the s-type and the d-type as well as the f-type pairing 
potential, and for the r^ parameter in the range of 1- 
10 a.u. Decreasing the density from rg — l to rs=3, all 
energies decrease by two orders of magnitude, from a 
few meV to a few per cent of meV. Further dilution of 
the electron gas to rs = 10 lowers the energies by another 
two orders of magnitude. This shows how delicate the 
balance is between the superconducting phase and the 
normal state. 

Slightly negative values of the condensation energies 
at r^^lO are obtained for f-waves and higher angular 
momentum pairing, but all components are very small 
and the most negative value is of order ~1.5T0~^ mcV. 
Since the biggest contributions to this negative energy are 



J con 



Tfc: 



5^ 0.0004 
fl 0.0002 




2 4 

angular momentum 



FIG. 2: The condensation energy, f condensation, and its com- 
ponents: the difference between the superconducting and the 
normal state of the RPA correlation energy /rpa (S-N), and 
that of the exchange energy fx (S-N), and of the kinetic en- 
ergy Tkin (S-N) and the anomalous Hartree energy /ah (this 
is nonzero only in the superconducting state). All energies are 
shown calculated at two densities: rs=l a.u. (upper panel) 
and and rs = 10 a.u. (lower panel). The model parameters in 
the pairing potential were fixed at 5 = 0.01- jj, and a — O.l-kp. 
Zero temperature was assumed. 



1 4.71e-0 6.56e-0 1.48e-0 7.95e-0 1.84e-0 

3 2.30e-2 9.37e-2 1.81e-2 9.79e-2 7.59e-2 

5 2.82e-3 1.28e-2 2.33e-3 1.27e-2 5.90e-4 

7 6.81e-4 3.33e-3 6.05e-4 3.29e-3 l.lOe-4 

8.5 2.87e-4 1.55e-3 2.78e-4 1.52e-3 4.16e-5 

10 1.19e-4 8.35e-4 1.45e-4 7.91e-4 1.84e-5 

d-waves 

1 8.28e-l 1.64e-0 3.93e-l 1.56e-0 4.43e-l 

3 2.04e-3 2.18e-2 4.80e-3 1.92e-2 1.82e-3 

5 4.21e-4 2.82e-3 6.19e-4 2.48e-3 1.42e-4 

7 8.40e-5 7.49e-4 1.61e-4 6.46e-4 2.64e-5 

8.5 2.16e-5 3.59e-4 7.38e-5 2.97e-4 9.99e-6 

10 7.80e-7 1.97e-4 3.85e-5 1.55e-4 4.43e-6 



f-waves 



1 7.24e-l l.lOe-0 3.02e-l 1.19e-0 3.32e-l 

3 3.Q7e-3 1.56e-2 3.69e-3 1.36e-2 1.37e-3 

5 3.31e-4 2.01e-3 4.75e-4 1.76e-3 1.06e-4 

7 5.87e-5 5.41e-4 1.23e-4 4.57e-4 1.97e-5 

8.5 1.21e-5 2.62e-4 5.66e-5 2.10e-4 7.48e-6 

10 -1.45e-6 1.44e-4 2.95e-5 l.lOe-4 3.32e-6 

TABLE I: The condensation energy, fcond., and its compo- 
nents: fRPA (S-N), /, (S-N), fAH and Tkin (S-N); calcu- 
lated for s- waves, d-waves, and f-waves for fixed 5— 0.01- jj, 
and a—O.l-kp, and at zero temperature (T=0 K) for differ- 
ent electronic densities, Vs- Energies are given in meV, and 
parameters: /i, kp and Vs are in atomic units. 



Jrpa and fAH, which are of order 1-10^* nieV while the 
condensation energy is of order 1-10~® meV, and since 
we trust to our numerical results up to the three leading 
digits, it is plausible that a change by one in the last 
position in the correlation and the anomalous Hartree 
energy may cause the change of a sign in the condensation 
energy. Actually for the same reason, the condensation 
energy at rs=10 for the d-wave pairing could be negative 
(because it is three orders of magnitude smaller than the 
leading contributions). But this uncertainty due to the 
numerical accuracy will happen neither for s-wave pairing 
nor for smaller r^ parameters, as one can see in Table 1. 
Later, we will show calculations for r^^lO with another 
choice of 5 and a. 

Finally, we change the parameters in the model pairing 
potential. These results are presented in Table 2 for Ts—I 
and for the angular momentum up to 3 (since before we 
found the possibility for a phase transition for f-waves). 

Firstly, we change a parameter for 0.05-fcF (before it 
has been fixed at Q.\-kp), and we keep the same 5 as in 



Table 1, i.e. 0.01-/J,. For this choice of 5 and /i, all ener- 
gies are smaller, and we do not find superconductivity at 
rs = l up to /— 3. The typical band widths at the Fermi 
surface are much more narrow than 0.1-kp. The results 
of calculations for s- waves performed by previous authors 
|25j | were smooth in a, and our results with two values of 
a appear smooth with respect to a variation of the angu- 
lar momentum. Therefore, we do not expect any change 
in conclusions by changing the a parameter. 

Secondly, we change 5 parameter for 0.001-^ (before it 
has been fixed at 0.01 -/x) and we keep the same a as in 
Table 1, i.e. O.l-fci^. Now, all energies are smaller, and 
the decrease of the condensation energy is about one and 
half order of magnitude while the change of the pairing 
amplitude (5 is by a factor of 0.1. The BCS behavior of 
the condensation energy is proportional to a square of the 
gap, ~ A^. We find in our calculations that, the depen- 
dence of the condensation energy on the gap amplitude 
is a bit weaker than the BCS one. For the lower momen- 
tum channels a power of that dependence is smaller than 



^ fcond. ^/r 



J- ki: 



rs = l, S=O.OhkF, (7=0. 05-^ 

3.22e-0 4.91e-0 1.09e-0 5.65e-0 1.39e-0 

1 9.18e-l 2.03e-0 4.47e-l 1.97e-0 5.37e-l 

2 5.39e-l 1.34e-0 3.05e-l 1.23e-0 3.52e-l 

3 3.94e-l l.OOe-0 2.37e-l 8.90e-l 2.66e-l 

rs = l, 5=0.001-kF, g-^O.l-M 

1.07e-l 2.24e-l 3.90e-2 2.58e-l 3.39e-2 

1 2.59e-2 8.12e-2 1.44e-2 8.06e-2 1.20e-2 

2 1.43e-2 4.97e-2 9.23e-3 4.72e-2 7.53e-3 

3 1.02e-2 3.53e-2 6.90e-3 3.31e-2 5.53e-3 



rs = 10, S=0.1-kF, a=0.5-fi 



7.85e-3 2.90e-2 5.05e-3 3.04e-2 1.40e-3 

1 1.49e-3 1.03e-2 2.23e-3 9.02e-3 5.39e-4 

2 1.13e-3 5.91e-3 1.58e-3 5.11e-3 3.53e-4 

3 9.17e-4 4.05e-3 1.25e-3 3.45e-3 2.67e-4 

TABLE II: The condensation energy and its components: 
/rpa (S-N), fa: (S-N), fAH and Tkin (S-N), calculated as func- 
tions of the angular momentum I for three sets of parameters: 
Vs, 5 and a. Energies are given in meV. 



for the higher angular momentum channels. This seem 
to be in contrast with measurements for ordinary, BCS 
like, superconductors which have a gap of the s-type. On 
the other hand, measurements are not able to split the 
purely electronic and the phononic contributions, and we 
did not add phonons to our calculations. 

Our result indicates that, wc should not expect super- 
conductivity at small r^ . Since all the energy components 
grow with 5 and a, we can say that the maximally neg- 
ative condensation energy would be for a=kp (all states 
contribute to pairing) and that there is no upper limit for 
S. This is due to the fact that our condensation energy 
is proportional to A" with 1.5<n<2. The value of 5, 
which we set to 0.01-^, is much bigger than typical gaps. 
For instance for Nb, we have s-wave pairing, rs=0.87, 
^=33.13 eV, while the experimental gap is 1.55 meV. We 
have chosen large delta, i.e. 0.01-/i~330 meV, for most of 
our calculations for the sake of accuracy, since the conclu- 
sions about the angular dependence of the condensation 
energy do not change with this parameter. 

For lower density, i.e. r^^lO, we calculated again the 
condensation energy but for a new choice of the gap pa- 
rameters: (5=0. 1-/^ and a^0.5-kp, which are unrealisti- 
cally high but the numerical accuracy is much better in 
this case. The result of above calculations is negative for 
superconductivity, which makes also a situation, that the 
phase transition occurs for a small gap, less probabilistic. 

If we wanted to include the electron-phonon or 
electron-paramagnon interaction, the expression for the 



total energy would be frequency independent because the 
frequency is integrated out in the SCDFT scheme (see for 
instance Ref. [isl lltil|V Similar to the total energy, also 
the gap function would be static, as it is now. This is in 
contrast to the Eliashberg formalism where inclusion of 
strong cou plin g changes the gap function to a dynamic 
parameter [30| ■ As for the feedback effect of phonons or 
spin-fluctuations for the electronic energy and vice versa, 
this effect would exist if one solved the SCDFT equations 
in the Bogoliubov-de Gennes forrn (for this formulation 
of SCDFT see for example Ref. [23 )• In the way how 
we calculate the condensation energy in this work, the 
aforementioned feedback would not exist since all energy 
components contribute to the total energy independently 
through a single equation. This is another difference from 
the Eliashberg scheme, where three coupled equations 
have to be solved for the self-energy S, pairing- vertex <& 
and polarization 11 |29l l30| . 

Summarizing results: the condensation energy is posi- 
tive for all angular momentums, which indicates that we 
are unlikely to obtain superconductivity from the elec- 
tronic correlations only. For the very dilute gas at rs = 10, 
we obtained slightly negative values of the condensation 
energy when the gap is small. These results should be 
however viewed with care, because all the energy com- 
ponents are very small and at the limit of numerical ac- 
curacy for the multidimensional quadrature. We trust to 
three leading digits in Tables 1 and 2, while the conden- 
sation energy at rs=10 is 3 orders of magnitude smaller 
than the biggest contributions. Also the random phase 
approximation d2j for the correlation energy is exact 
only in the limit of high densitv. pOl| i.e. rs<l. This treat- 
ment of the Coulomb interaction might be insufficient for 
the density at rs = 10. Whether it is an appropriate ap- 
proach it depends also on the angular momentum of the 
pairing potential and on the physical property one is in- 
terested 111 mi. 



IV. DISCUSSION AND CONCLUSIONS 

Whether the superconductivity can exist without 
phonons or not is a very old problem. In 1965, Kohn 
and Luttinger p suggested a mechanism of the Cooper 
pair a\ formation in the homogeneous gas due to Friedel 
oscillations 7] . These authors did not assume any partic- 
ular form of the interaction, which could even be purely 
repulsive, since the attractive regions could form in real 
space because of a sharpness of the Fermi surface in the 
reciprocal space. It has been discussed that, for the pair- 
ing potential at higher angular momentum, the super- 
conducting state was more favorable than for s- waves. 
The above conclusions were based on the mathematical 
analysis of irreducible vertexes with the particle-particle 
interaction up to the second order (see Fig. QJ . The cri- 
terion used for the superconductivity was the occurrence 
of a pole at Tc in the scattering amplitude for pairs of 
quasi-particles of equal and opposite momenta and in the 



total energy corresponding to two particles on the Fermi 
surface. 

Later work on the superconducting homogeneous gas 
within the Eliashberg theory, treating the Coulomb inter- 
actions on the RPA level and beyond, can be compared 
to the results presented in this paper. 

Superconductivity obtained due to the plasmon ex- 
change alone seemed to be overestimated. Several papers 
solving the Eliashberg equations with the RPA model 
for the Coulomb interaction predicted superconductiv- 
ity for s-wave pairing at quite high densities. Within 
the weak coupling limit of the electron-phonon inter- 
action the Eliashberg equations could be considered as 
k-dependent only due to the KMK approximation y\. 
In the limit of strong coupling to phonons, the k- and 
cj-dependent equations have to be solved. While first 
approximation led to superconductivity .q at the den- 
sity lower than rs=6, solving the strong coupling regime 
equations yields a change in a sign of the Coulomb pa- 
rameter (At*<0) at the density rs~2.5. The total 
energy diagrams included in both approaches are shown 
in Fig.n The common assumption in the aforementioned 
two approaches is that the polarization function, which 
enters the RPA screened interaction contains only normal 
Green's function loops, neglecting the anomalous ones 
which are also included in our scheme. In addition, the 
normal Green's function used by Sham and co-workers 
was obtained self-consistently in contrast to all the other 
papers discussed here. The correctness of such an ap- 
proach was discussed by several autho rs Jlo . 14^ also by 
occasion of the GW approximation *43!j44|. The criterion 
for superconductivity used by Takada g] was a nonzero 
critical temperature. Sham and co-workers [lOl \13L IliL l4a | 
considered the electron gas to be superconducting when 
the Coulomb pseudo-potential fi* was positive. 

A number of papers that included vertex corrections on 
different level approximations predicted decreased tem- 
peratures of phase transition. Grabowski and Sham [l3| 
studied the w-dependent only Eliashberg equations with 
vertex corrections up to the second order. Within that 
simplified model, the signum of /i* changed at rs=7. 
More extensive study of vertex corrections were done 
within the KMK, k-dependent scheme extended to the 
strongly correlated systems by Takada [I^ , who included 
more than 50 diagrams systematically using the effective- 
potential expansion. As a conclusion of that work the 
phase transition occurred at rs>3.9, the maximum of 
the critical temperature was obtained at rs=7.2, an d Tr 
decreased for lower densities. Biiche and Rietschel [lj| 
added the vertex corrections, within the phenomenologi- 
cal model by Kukkonen and Overhauser (KO) 45], to the 
earlier work by Rietschel and Sham 10], and did not ob- 
tain superconductivity up to rs=5. The parameter fi* for 
gas in the range of l<rs<5 was positive and varied be- 
tween 0.05 and 0.1. Then, Takada [T^ performed k- and 
w-dependent calculations with the local-field correction 
of the KO model and showed a significant effect of cor- 
rections beyond RPA around rs=5. Another observation 



made in his work was that, although the compressibility 
K and the spin susceptibility x were strongly dependent 
on the vertex corrections, Tc of gas at the density rs>20 
was similar to the temperature obtained from RPA. For 
rs>40 the critical temperature has been approximated 
by Tc^OM-Ep. 

All the aforementioned papers on the superconducting 
homogeneous gas, within RPA and beyond, dealt with 
the pairing potential of the s-type. The s-, p- and d-type 
pairing potentials examined in Ref. |4l| led to the con- 
clusion that the vertex corrections to s-waves are much 
more important than to the higher momentum channels. 

In this work, we calculated the condensation energy at 
zero temperature, instead of solving the gap equation for 
finite temperatures. This way we are not able to deter- 
mine the critical temperature for the densities at which 
we find superconductivity. On the other hand, negative 
values of the condensation energy are very small, of order 
^1.5T0~^ meV, at the densities rs>10 for f-waves and 
the higher-? channels. For d- waves, the condensation en- 
ergy is positive but very small. 

We concluded that this effect might be due to an ac- 
curacy of the multidimensional numerical quadrature of 
objects with many singularities. On this point we want 
to comment that most of standard quantum chemistry 
programs which calculate two-electron integrals do not 
accede an accuracy higher than 6 important digits. Our 
task is even more difficult because in addition to cal- 
culating the k- and q-momentum vectors, we have to 
perform the quadrature over the bosonic Matsubara fre- 
quency, and singularities in the superconducting state 
have much more complicated shape |3a ] than those of the 
two-electron integrals calculated for the normal state. 

The positive aspect of the method employed here is 
the absence of any approximation except RPA. The DFT 
is exact for the ground state studied here. We do not 
make any assumption about phononic interactions, which 
we neglect. But if we wanted to include the electron- 
phonon interactions, then the way of treating the weak 
and the strong coupling would be the same |22]. We 
do not drop either the momentum dependence or the 
frequency dependence of the Coulomb interaction. The 
vertex corrections, especially important for the moderate 
densities for the properties like the critical temperature 
[iq]. are not taken into account in this work. However, we 
believe that the inclusion of vertex corrections would not 
change the conclusions, because for f-waves the strength 
of the Coulomb attraction is not as much overestimated 
by the polaron exchange as for s-waves [IJ . 

Another interesting question would be whether it is 
possible that at some densities s-waves are favorable for 
the superconducting state and at other densities the 
higher-Z pairing would lead to lower energy. Such s- 
wave to p-wave transition has been reported by Takada 
|4ll ] at rs=4.7 while for higher densities the energy of p- 
waves was lower. KiichenhofF and Wolfle |4a], by solving 
two coupled Bethe-Salpeter equations for the two-particle 
vertex functions in the particle- hole channels, found the 
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p-wave superconductivity for 10<rs<35 and the s-wave 
superconductivity for rs> 35. From our results, which are 
monotonic with the angular momentum number, such l- 
wave to Z'-wave transition seems not to be the case [M |. 

We assumed singlet pairing. Thus, only the even num- 
bers I (like s-waves, d-waves etc.) are relevant, and it 
would be incorrect from the symmetry point of view to 
compare the total energy of s-waves with the total en- 
ergy of p-waves. On the other hand, if there is not 
much energy gain by the Cooper pair formation then the 
spin pairing in the superconducting phase probably is the 
same as the magnetic phase of the normal state. Within 
the QMC methods J43| , it has been widely examined the- 
oretically that the ground state of the homogeneous elec- 
tron gas is paramagnetic for high and intermediate densi- 
ties, and the transition to the ferromagnetic phase occurs 
at about rs~25. There is an experimental evidence |43| 
for the ferromagnetic phase in Cai_a;Laa;B6 at the den- 
sity rs=28 a.u., where the saturation moment of 0.07 fiB 
per electron resists below the temperature 600 K, which 
is of the order of the Fermi temperature of the electron 
gas. However, for this experiment the iron substrate was 
chosen. Another experimental group |49j| could not find 
any evidence for the intrinsic ferromagnetism in any of 
Ai_a;Laa;Be (A=Ca, Sr) samples. As for the novel su- 
perconducting materials with triplet pairing, they can- 
not be described within a model based on the local spin 
density approximation, where the parametrization on a 
homogeneous gas results works well. In order to describe 
those superconductors, like Sr2Ru04 or (TMTSF)2X, it 
has been shown by Shimahara in Ref. [5Q | that, one 
should add strong short range correlations to weaken the 
Coulomb interaction and then the electron-phonon mech- 
anism could cause pairing. Moreover a model assumed 
in our work is 3D, while Sr2Ru04 should be described in 
2D and (TMTSFJzX in ID. 

There is also a question about an effect of the fluctu- 
ations which could mediate the pairing interaction (for 
such a model see the work of Abanov and Chubukov in 
Ref. 51] ). Such calculations for the condensation en- 



ergy in strongly correlated systems, where both effects of 
the electron-phonon and the electron-paramagnon inter- 
actions are taken into account, have been performed by 
Haslinger and Chubukov in Ref. 113 • The essential differ- 
ence between those calculations and ours, if we included 
phonons and paramagnons, would be in the frequency de- 
pendence of the gap function, which in our case is static 
due to different formulation of the problem from the very 
begin. For the inclusion of the dynamical effect to the 
gap, one needs to go beyond the Born-Oppenheimer ap- 
proximation. In fact, such formulation exists within the 
multicomponent DFT scheme proposed by Kreibich and 
Gross [5^, and developed by Van Leeuwen [53|. 

In conclusion: we have calculated the condensation en- 
ergy of the homogeneous electron gas at zero temperature 
within the density functional theory for superconductors. 
The random phase approximation was assumed for the 
Coulomb interaction and no phononic contributions have 
been added. We did not consider pairing mechanism me- 
diated by fluctuations. Within this approach, there is 
no superconductivity for any momentum of the pairing 
potential for the densities up to rs~9. We found very 
weak superconductivity for f-waves and higher-/ pairing 
at Ts—lO, but this effect is so small that could be due 
to the neglecting of the vertex correction, or due to the 
accuracy of the numerical quadrature. 
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